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Abstract — This paper considers the problem of time-domain 
assessment of the Phase Margin (PM) of a Single Input 
Single Output (SISO) Linear Time-Invariant (LTI) system using 
a singular perturbation approach, where a SISO LTI fast 
loop system, whose phase lag increases monotonically with 
frequency, is introduced into the loop as a singular perturbation 
with a singular perturbation (time-scale separation) parameter 
e. First, a bijective relationship between the Singular Pertur- 
bation Margin (SPM) e max and the PM of the nominal (slow) 
system is established with an approximation error on the order 
of £ 2 . In proving this result, relationships between the singular 
perturbation parameter e, PM of the perturbed system, PM 
and SPM of the nominal system, and the (monotonically 
increasing) phase of the fast system are also revealed. These 
results make it possible to assess the PM of the nominal 
system in the time-domain for SISO LTI systems using the 
SPM with a standardized testing system called ”PM-gauge,” as 
demonstrated by examples. PM is a widely used stability margin 
for LTI control system design and certification. Unfortunately, it 
is not applicable to Linear Time- Varying (LTV) and Nonlinear 
Time- Varying (NLTV) systems. The approach developed here 
can be used to establish a theoretical as well as practical 
metric of stability margin for LTV and NLTV systems using a 
standardized SPM that is backward compatible with PM. 


I. INTRODUCTION 

Due to complexity of the dynamics of practical systems, 
feedback controllers are designed almost invariably based 
on simplified but concrete mathematical models known 
as nominal systems under some simplifying assumptions. 
However, the designed controllers will be implemented on 
the actual systems, which differ from the mathematical 
nominal model in many uncertain ways that are known as 
(structural) perturbations to the nominal system. The actual 
system or any mathematical model with higher fidelity than 
the nominal model are then called the perturbed system. In 
nonlinear control theory, perturbations that do not change the 
order of the nominal system, such as parametric dispersions 
or time-dependent variations, and negligible nonlinearities, 
are called regular perturbations; whereas those that change 
(increase) the order of the nominal system, such as fast 
parasitic dynamics that are ignored during controller design, 
are known as singular perturbations. An effective controller 
must be able to guarantee stability of the actual system 
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in the presence of these perturbations. References [l]-[3] 
present an overview of singular perturbations including 
perturbation models, numerical and analytical treatments, 
ideas, results, and the applications of singular perturbation 
theory in different fields. 

The extent to which a controller accommodate structural 
perturbations is called stability margin. Note that stability 
margin is different from the concept of domain of 
attraction, which is the ability of the controller to cope 
with perturbations in initial states. It is also different from 
the concept of input-output sensitivity (input-output gain), 
which indicates the ability of the closed-loop system to 
attenuates perturbations in the input, such as disturbances, 
infeasible tracking command trajectories or noises. Based 
upon the perturbation categories, stability margin can be 
divided into singular and regular perturbation margins. For 
Single-Input-Single-Output (SISO) Linear Time-Invariant 
(LTI) systems. Gain Margin (GM) and Phase Margin (PM) 
have been successfully used to gauge a nominal system’s 
capability in accommodating parametric uncertainties 
(regular perturbation) and parasitic dynamics (singular 
perturbation). While parametric perturbations do affect the 
PM, PM is predominantly a margin for accommodating 
parasitic dynamics as a singular perturbation margin. 
Whereas parasitic dynamics (with non-unity gain) do 
affect the loop gain, GM can be considered predominantly 
a regular perturbation margin. References [4] -[6] have 
discussed Nyquist robust stability margin, time-delay 
stability margin and so on, in the fields of robust control 
and L\ adaptive control including regular and singular 
perturbation margins from the view of above classification. 

However, since PM and GM are based on frequency 
response of LTI systems, they are not applicable to Linear 
Time-Varying (LTV), Non-Linear (NL) systems (including 
NLTI and NLTV) systems, which do not lend themselves to 
frequency domain analysis or synthesis. A common practice 
in controller design for LTV and NL systems applies the 
LTI controller design methods at discrete operating states 
within the operating envelope or along a nominal trajectory 
under the assumptions that the nonlinearity is benign 
and time-variance of the system parameters or nominal 
trajectory are slow. In operation, the controller gains are 
’’scheduled” according to the current operating state or 
commanded trajectory. In essence, such gain-scheduled 
controllers and the closed-loop system are NLTV in nature. 
But the current practice is to certify such controllers by the 



GM and PM at the design points, which is well-known to 
be neither necessary nor sufficient for stability should the 
assumptions of benign nonlinearity and slowly-varying fail 
in operation, such as departure or upset, or simply flying a 
maneuvering trajectory such as in a flight control system, 
thereby jeopardizing the safety of the system. The benign 
nonlinearity and slowly-varying assumptions also limit the 
performance of the control system. 

In recent years many advanced nonlinear and time-varying 
control techniques have been developed for NLTV systems 
to alleviate the benign nonlinearity and slowly-varying 
assumptions, thereby fully exploiting their performance 
potentials and improving their safety. However, there does 
not appear to have been reported for NLTV systems any 
stability margin metrics that are: (i) theoretically based, 
(ii) practically measurable, (iii) backward compatible in 
the sense that when applied to LTI systems, these metrics 
have a bijective correspondence with the GM and PM [7]. 
The lack of such effective stability metrics has inevitably 
hampered the acceptance of advanced LTV and NL control 
techniques in control engineering practices. 

In this paper, we propose a time-domain approach for 
assessing PM of SISO LTI systems based on the singular 
perturbation theory, which can potentially be extended to 
LTV and NL systems as an effective gauge of stability 
margin for parasitic dynamics similar to what PM does for 
LTI systems. 

The research on singularly perturbed systems with two- 
or multi-time scale are of theoretical and practical interests 
[8]- [10]. When the general form of the LTI autonomous sin- 
gularly perturbed system is considered, it can be represented 
as [11] 

x(t) = A\\x(t) +Ai 2 z(t) 
ez(t) = A 21 x(t)+A 2 2 Z(t) (1) 

where x(t) £ R",z(t) £ R m are state vectors of slow and 

fast dynamics (two-time scale), respectively; An, A \ 2 , A 2 \, 
A 22 are constant matrices and A 22 is Hurwitz; e is a small 
positive parameter that describes the time-scale separation, 
such as the ratio between the bandwidths or between the 
time-constants of the slow and fast dynamics. When £ = 0, 
(1) reduces to the nominal system 

x(t) = [An — Ai 2 A 22 A 2 i]x(t) (2) 

In [12] and [13], the perturbed system (1) was decomposed 
with Chang transformation [14] into two decoupled 

subsystems of the slow and fast dynamics. This result 
was generalized in [15] (see also [16]) to multi-time-scale 
singularly perturbed linear systems. The multi-time-scale 
decomposition can be used for implementing linear dynamic 
controllers, observers, and Kalman filters using different 
sampling rates. 


Many results on stability of singularly perturbed NLTV 
systems are available (cf. Chapter 11 of [11] and Tikhonov’s 
Theorem, p434 of [11]), which, roughly speaking, establish 
that when the null (quasi) equilibrium states of both the 
nominal operating error dynamics (defined by the slow 
system trajectory relative to the nominal operating state or 
trajectory) and the boundary layer dynamics (defined by 
the fast system trajectory relative to the slow dynamics) 
are exponentially stable, there exists a constant e max > 0 
such that the null equilibrium of the singularly perturbed 
system (the composite system) is exponentially stable for all 
e < e max . We may call e max Singular Perturbation Margin 
(SPM). While both SPM and PM provide an indication 
of the stability margin of the nominal (slow) system for 
parasitic dynamics, the quantitative relationship between 
the SPM e max and the PM does not appear to have been 
studied. 

The main results of this paper offer a singular perturbation 
approach for time-domain assessment of phase margin. Here, 
the nominal system (2) is assumed to have been (exponen- 
tially) stabilized with a healthy PM, and the fast system can 
be written as 

£z(t) = A 22 z(t) + Bfastx(f) 

y (0 = Cfastz(0 (3) 

where Cf as t = A\ 2 and Bf ast = Agi. Fast system (3) is 
assumed to be a monotonic-lag system (i.e. phase lag 
increases monotonically with frequency) with sufficiently 
large phase lag (larger than the PM of the nominal system) 
at the gain cutoff frequency ft) C g,nom of the nominal system 
for sufficiently small e. It will be shown that the PM of the 
perturbed system can be approximated by the sum of the 
PM of the nominal system and the phase of the fast system 
at ft) C g,nom of the nominal system with an error on the order 
of £ 2 . With the ft) C g,nom and £ max , our Theorem 1 establishes 
a bijective relationship between the SPM £ max and the PM 
of the nominal system with an approximation error on the 
order of £ 2 . This result can be used to assess the PM of the 
nominal system in the time domain when the SPM £ max is 
found either by direct calculation of the eigenvalues of the 
perturbed system, or estimated using a singular perturbation 
approach. In practice, the fast system is designed as a ”PM 
gauge,” which should be optimized and standardized in 
order to obtain a quantitative estimate of PM of the nominal 
system. Proof of Theorem 1 revealed that the PM estimation 
error is due to the magnitude frequency response of the fast 
system at ft) C g.nom, therefore the optimal PM gauging system 
should be an allpass filter with unity gain for all £ < £ max . 
The usage and effectiveness of this method are shown by 
examples. 

The significance of the main results of this paper is not 
merely an alternative method for determining the PM of a 
SISO LTI system; rather it provides an avenue for defining a 
stability margin for LTV and NL systems that is theoretically 



based, practically measurable and backward compatible 
with PM. The extension relies on theoretically rigorous and 
practically computable algorithms for computing the SPM 
for LTV and NL systems. For LTV systems, a recently 
developed differential algebraic spectral theory [17]-[22] 
could be used to assess the SPM in a manner similar to the 
aforementioned Chang transformation and multi-time-scale 
decomposition for LTI systems. The LTV SPM methods, 
once developed, could be extended to NLTV systems by 
virtue of the equivalence of exponential stability of an 
equilibrium state of a NL system and its linearization 
thereof. For NL systems (including LTV systems as 
a special class), well-known Lyapunov function based 
singular perturbation results (see, for instance. Theorems 
11.3 and 11.4 in [11]) can be used to obtain conservative 
estimates of SPM. In this latter approach, a recent result 
on spectral Lyapunov function [23] would prove to be useful. 


and e max is defined to the smallest ft)BW,fast that makes the 
perturbed loop system be the £ at onset of instability. The 
magnitude M(a>) of transfer function L is defined by 

M(co) = 20 log \L(jco)\ 

The perturbed system (1) is shown in Fig.l. 



Fast System 


The paper is organized as follows. After this introductory 
section, our main theoretical results will be presented 
in Section 2. The singular perturbation approach for 
time-domain assessment of PM is illustrated by examples 
in Section 3. Section 4 concludes the paper with some 
insightful remarks and suggestions for further research. 


II. MAIN RESULT 

The symbols used in the theorem statement and proof are 
shown in TABLE I. 

TABLE I 
Symbols 


Symbols 

Meanings 

Unit 

PM pt 

Phase margin of the perturbed system 

rad 

P Mn om 

Phase margin of the nominal system 

rad 

PM fast 

Phase margin of the fast system 

rad 

0pt(®,e) 

Phase of the perturbed system at freq. co, 
with singular perturbation parameter e 

rad 

0 nom(ffl) 

Phase of the nominal system at freq. co 

rad 

0fa s t(®,£) 

Phase of the fast system at frequency co, 
with singular perturbation parameter e 

rad 

M pt (co,e) 

Magnitude of the perturbed system at freq. 
co, with singular perturbation parameter £ 

dB 

Miom (&)) 

Magnitude of nominal system at freq. co 

dB 

Mfast(ft>,£) 

Magnitude of the fast system at freq. co, 
with singular perturbation parameter £ 

dB 

fflcg.pt 

Perturbed system gain crossover frequency 

rad/sec 

^cg,nom 

Nominal system gain crossover frequency 

rad/sec 

^BW,nom 

Bandwidth of the nominal system 

rad/sec 

®BW,fast 

Bandwidth of the fast system 

rad/sec 

^pt (^5 £) 

Perturbed system open loop transfer function 
with singular perturbation parameter £ 


^fast(^? 

Fast system open loop transfer function, 
with singular perturbation parameter £ 



The singular perturbation parameter £ is defined by 
tAiW'.nom fr) C p,nom 

£ = S3 

&*BW,fast t^BW.fast 


Fig. 1. Singular Perturbed System 

The transfer function of the fast system (3) is given by 
Cf a st(v = Cf as t [£.?/ A 22 ] Bf as t (5) 

A. Theorem 1 

Consider the singularly perturbed LTI system (1), where 
x{t) £ R", z(t ) £ R m ; An, A 12 , A 21 , A 22 are constant matrices 
and A 22 is Hurwitz; £ > 0 is the perturbation parameter. If 
the fast system is a SISO lag system with unity DC gain and 
monotonic phase lag that is sufficiently large for sufficiently 
small £, then PM nom and £ max satisfies a bijective relationship 
defined by 

P Aurora = 0lasl ( PLg.norn - Fmax ) | V 0(£~) (6) 

Proof of Theorem 1 relies on the following lemma, which 
is an important result in its own rights. 

B. Lemma 1 

Consider the singular perturbed LTI system (1), with 
the same assumptions and definitions of the parameters in 
Theorem 1. The following relation between PM pt , PM nom 
and (j) f as t is satisfied 

PM p t — PM nom 0f a st ( Ptcg.nom ; £) T D ( £ ) (7) 


C. Proof of Lemma 1 

Suppose A 22 is Hurwitz, the loop transfer function of the 
fast system, which is a rational function with unity DC gain, 
and poles having negative real parts, can be written as 

L {o) ngi [ej/zi + 1] nf 2 | [(esM„,) 2 + (2C,„,-M„,)£S + 1] 

U I jL l [£s/pj+l}U%i[(£s/(0, lj ) 2 + (2^n j /c0n j )es+l\ 

The fast system is a phase lag system, so 
min { {z,}" 1 ! U { CO IHl } > min { {pj}% U { } (8) 


and 


Let PMpt — PM nom 4" 0fast(&)cg,nom)> ^.nd 


< 


M\ M2 

ni es Ai + linN-K) 2 + (2 ^ ni /(Om,)es+\\ 

i= 1 i'=l 

JVi IV 2 

n^M- + 1] nt( es /^) 2 + (2^ nj /co nj )Es+ 1] 


;=1 


7=1 


Let 


(Bmin = mill {{pj}j = \ U } 

and for the magnitude of Lf ast , we have 

|^Kfast(®cg ,nom) | 


< N\ 


20 log 


; Klcg.nom , , 

J Wmin + 


+N 2 


20 log 


1 + Ej- 


ZWcg ,nom 1 ( q ' ^cg.nom \ 2 
^min ®min 


(9) 


where £ £ (0,1). By Taylor Expansion of (9) at the point 
e = 0 , we obtain 


\PM pt -PM pt \ 

= | (k + 0nom (t^cg,pt) H“ 0fast(^cg,nom)) 

(P M nom + 0fast(tt)cg,nom)) | 

= | (?£ 4" 0nom (®cg,pt) + 0fast(ft>cg,nom)) 

( ^ 4" 0nom (G7 C g ,nom ) H“ 0fast(t^cg,nom)) | 

^ | 0nom ( t^cg,pt ) 0nom ( &)cg,nom ) | 

4" | 0fast(^cg,pt) — 0fast(t^cg ,nom ) | 

= |Aft)| ( 4_(|)nom(©l) + ^fe St (® 2 ) ) (14) 

for some ( 0 \ , CO 2 G [min{co C g 5 pt, &)cg,nom} , max-[n)cg,pt 5 &)cg,nom}]> 
and then 


\PMpt~ PM pt \ = 0(|A«b|) = 0(e 2 ) 


(15) 


Therefore 


PM- pt — PM n oxn T 0fast(tt)cg,nom) T 0{£ ) 


|^fast(^cg ,nom) | 
2 


< 2(Wi 


-ft), 


cg,nom 


2(11110)0), 


2 1 

min 


0 (e 4 ) 


+2(W 2 

= 0(e 2 ) 


(1 2£ )ft> cg , nom 2 

0"io)<C e “ 


0 (e 4 ) 


(10) 


Since M nom (ft) C g, n0 m) = 0, the magnitude of the perturbed 
system at ft) cg , n om is 


T/pi ( ft)cg ,nom ) — M nom ( C0 C g ,nom ) 4' Mf as t((o C g ,nom) 

= ^ffast ( tWcg,nom ) 

= 0 {£ 2 ) ( 11 ) 


which is shown by Fig. 2. 



Fig. 2. Magnitude-Frequency Plot 


Let | Aft) | — | C 0 C g ? pt ft) C g 9 nom 1 9 
theorem, we have 

|-^fpt(tWcg,nom) | 

| Aft) | 


and by the mean value 




( 12 ) 


for some coq G [niin^cOcgjpt^^cg^oml^Pdcix-^uJcgpt^cOcg^om}]* It 
then follows from ( 12 ) and the last equality in ( 11 ) that 




l^fast (t^cg ,nom) | 
|^M pt (fi)o)| 


0 ( 8 2 ) 


03) 


D. Proof of Theorem 1 

By assumption, the fast system has a monotonically in- 
creasing phase lag. When PM nom is given, |0f ast (ft)cg,nom,£)| 
increases as e increases, and PM pt decreases. When PM pt 
reaches 0, we obtain £ max . By Lemma 1 

PM pt 0 PM nom + 0fast ( ft)cg,nom - £max) T 19 ( £ ) 

Therefore 

PM nom — 1 0fast ( ft)cg,nom , £max ) | T D(£“) 

Notice that 0f as t(s,£) = ^fastAs). It follows the assumption 
that 0 fast (s,e) is a monotonically increasing function of ft) 
that PM nom is a bijective function of the SPM £ max . 


III. Time-Domain Phase Margin Assessment 

In this section we will apply Theorem 1 and Lemma 1 to 
time-domain assessment of PM nom for the nominal system. 
We will assume the gain-crossover frequency ft) C g,nom is 
known for the given nominal system, which will be used to 
design a LTI singular perturbation (fast) system Lf ast (.v, £), 
called a PM-gauge, whose phase 0fastC/ft)cg,nom,£max) will 
be used to determine PM nom according to Theorem 1, 
where the SPM £ max can be determined in the frequency 
domain from the condition PM pt (£ max ) = 0. However, the 
purpose here is to calculate £ max in the time-domain from 
the state equations of the perturbed system comprising 
of the nominal and the PM-gauge systems. In Subsection 
3.1 we will first present both analytical and numerical 
methods for evaluating the SPM £ max in the time-domain. 
Then in subsequent subsections we will present three types 
of PM-gauge, each of which has unique advantages and 


disadvantages that are suited for a certain applications. 


The design methods and the effectiveness of the designed 
PM-gauges will be demonstrated using the nominal system 
shown in Fig. 3 with a loop transfer function given by 


i( s ) 


CpBpjKj + Kps) 
s(s + A P ) 


(16) 


which is a first-order LTI plant with a Proportional-Integral 
(PI) feedback controller, where the plant parameters are given 
by Ap = 1, Bp = 1, Cp = 1. In order to illustrate the 0(e 2 ) 
accuracy of estimated PM nom using different PM-gauges, two 
sets of controller gains K\ and Kp will be used, which result 
in two sets of PM nom and ft) C g,nom values, as shown in TABLE 
II. 


PI Controller PM-Gauge Plant 



where N = n + m. For e = £ max , (18) has an imaginary 
eigenvalue A = j/3, which simultaneously satisfies the real- 
part and imaginary-part of (18) 


( — l) <L fl2i+i/3“ <: + ... + fl5(e)/3 4 — c/3(e)/3“ + ai(e) — 0 

(-l^^/t+i)/^ + ■•■ + fl 6( £ )/3 4 ~ fl 4(£)/l“ + 02(e) = 0 

(19) 


where k = 1,2, ...,M and M = [A/2], and ajv+i(e) = 1, 
a,v 12 (e) = 0 when appropriate. Thus, a standard numerical 
method can be used to solve for /3 from both equations 
in (19) for a given e, and e = e max is determined by 
increment e until the values for ji from both equations 
are sufficiently close. To speed up the search process, 
increment of e can be variable based on the difference 
between the two solutions for /3 from both equations in (19). 

When the order N < 9, it is possible to solve for e max 
from (19) analytically using the root formulas for quadratic, 
cubic and quartic polynomials, as illustrated below in the 
PM-gauge design examples. 


Fig. 3. The Nominal System 


B. Minimum-Phase PM-Gauges 


TABLE II 

Nominal System Parameters 


Parameters 

Kp 

Ki 

P ^nom 

^cg.nom 

Case 1 

1 

3 

k/3 

1.7321 

Case 2 

0.5 

7.83 

it/ 6 

2.7320 


A. Time-Domain Calculation of SPM 

Consider the singularly perturbed LTI system (1) where 
A 22 is Hurwitz, and assume the nominal system (2) 
has been exponentially stabilized, i.e. [A n — Ai 2 Ai n 1 A 2 i] 
is Hurwitz. Then standard singular perturbation theory 
provides exponential stability of the perturbed system 
for sufficiently small e > 0. Thus, SPM e max can be 
defined as the least upper bound on e > 0 for which 
the perturbed system remains exponentially stable. The 
SPM e max can then be found by monotonically increasing 
e from a very small value until one of the eigenvalues 
of the perturbed system becomes imaginary (including zero). 


Because minimum-phase zeros introduce monotonic phase 
leads, which violate the monotonic phase lag assumption for 
Theorem 1, a minimum-phase PM-gauge must be an all -pole 
system. Since the magnitude function of an all-pole system 
is not a constant, a minimum-phase PM-gauge will have 
an 0(e 2 ) error in estimating PM nom according to Theorem 
1. The following example illustrates a minimum-phase 
PM-gauge. 


Consider a second-order minimum-phase PM-gauge with 
a transfer function 


^fast(^5 £) 


< 

(es) 2 + 2^coo(es) + ®o 


The phase function of Lf ast (yft), e) is given by 


|<fest(Me)| = arctan 



( 20 ) 


( 21 ) 


To this end, first rewrite the perturbed system (1) as a 
composite system 


x(t) 


z(0. 

1 

1 


-A21 

L£ 



x{t) 

_z(t) 


= A(e) 


x{t) 

z(0. 


The characteristic equation of (17) can be written as 


(17) 


det[A/— A(e)] = X N +aN(e)f N 1 + ... + 02 (e) h +a i(e) = 0 

( 18 ) 


We choose (0 0 = ffl C g,nom, so that e max = 1 corresponds to 
PMnom = tt/ 2 rad (90°). By choosing £ = 1 we have 


PM m 


l — |0fast(./^cg ,nom > £max ) | = arctari (7^H (22) 


with an 0(e 2 ) error. The approximate relationship between 
PM and SPM is plotted in Fig. 4. 



Fig. 4. Approximate Relationship between PM and SPM for (20) 


For the two test cases of nominal systems specified in 
TABLE II, the transfer function of the corresponding PM- 
gauges along with the corresponding SPM e max values, esti- 
mated PM nom values and the estimation errors, and the bound 
k= (PM n0 m Error)/e^ ax are shown in TABLE 2. The 0(£ 2 ) 
accuracy of the estimation error is clearly demonstrated. It is 
reassuring to know that the estimation error is significantly 
smaller as the design becomes more critical. The estimation 
error can be further reduced by using a higher order PM- 
gauge, such as cascading two PM-gauges of (20). However, it 
will only reduce the constant k in the last column of TABLE 
III, not by orders-of-magnitude. 

TABLE III 

Minimum-Phase PM-Gauge Test Results 


magnitude function of the PM-gauge Mf ast (s,£). Thus it 
is possible to design zero-error PM-gauges using all-pass 
filter transfer functions Lf ast (s,e) whose magnitude function 
Mf ast (s,e) = 1. To illustrate the design, consider first-order 
all-pass filter transfer function 

j £ S _ 

L {ast (s,e) = —§- (24) 

+ mo 


whose phase function is given by 


|<fest(Me)| = arctan 



(25) 


By choosing (Oq = fi) C g,nom, so that e max = 1 corresponds to 
PM aom =n / 2 rad (90°), we have 

PT/ no , n = 0| a st(./tO C g ,nom •> £max )l= ar ct an (^-) (26) 

' Lmx ' 

which is identical to (22). However, there is no estimation 
error, and the order of the PM-gauge is only half of that 
of the minimum-phase counterpart (20). The drawback of 
such zero-error PM-gauges is that they are necessarily non- 
minimum-phase, which may cause problem in simulation 
verifications of PM or SPM using such PM-gauges. The 
results of applying the zero-error PM-gauge (24) to the two 
test cases are documented in TABLE IV. 


TABLE IV 

Zero-Error PM-Gauge Test Results 


Parameters 


SPMe max 

Case 1 

3.0002 

(ef) 2 + 3 .4642 (e j) + 3 .0002 

0.7087 

Case 2 

7.4638 

(es) 2 + 5.4640(es) + 7.4638 

0.2819 


Param. 

Est. PMnom 

Exact PMnom 

P M n0 m Err. 

k 

Case 1 

1.2331 

~k]T==TX)4T2~ 

17.75% 

0.3701 

Case 2 

0.5495 

n/ 6 = 0.5236 

4.95% 

0.3259 


It is noted that the SPM e max determined in the frequency 
domain from PM pt (£ max ) and that determined using the 
time-domain method introduced above are identical. Since 
the order of the perturbed system is N = 4, e max can be 
determined analytically using the quadratic root formula 
from (19), which now becomes 

/3 4 -a 3 (e)/3 2 + ai(e) = 0 
a 4 (e)j3 2 + a 2 (e) = 0 

Eliminating /3 2 from these equations yields an equation for 
e 

a\{e) + a^{e)a 2 ,{e)a 2 {£) + a] t (e)ai{e) = 0 (23) 

which turns out to be a second-order polynomial in e 2 , and 
therefore can be solved using the quadratic root formula. 
Among the four roots, the smallest positive root is £ max . 

C. Zero-Error PM-Gauges 

From the proofs of Theorem 1 and Lemma 1 it is seen 
that the Off 2 ) estimation error is due to the non-unity 


Parameters 

Pfast( s i £) 

SPMe max 

Case 1 

1 

1 1.7321 

1 _i £L_ 

1 ^ 1.7321 

0.5774 

Case 2 

1 

1 2.7320 

1 _i £i_ 

1 ^ 2.7320 

0.2679 


Parameters 

Est. PMnom 

Exact PMnom 

P 4^ n0 m Err. 

Case 1 

1.0472 

7t/3 = 1.0472 

0 

Case 2 

0.5236 

n/6 = 0.5236 

0 


D. Approximately-Linear-Phase PM-Gauges 

Although the zero-error PM-gauge of the preceding sec- 
tion gives exact PM nom from SPM £ max , the relationship 
between them is a nonlinear function. In practice, it is 
desirable to have a linear relation between PM nom and SPM 
£ max with zero-error. A pure transport delay possesses the 
desired unity magnitude function and linear phase function. 
Thus, we may design linear-phase, zero-error PM-gauge as 

L fast ( S ,£)=e~ T » es (27) 

whose phase function is given by 

|#festC/ffl)£)| = Toes (28) 

By choosing 7o = 7t/0J C g J1()m , so that £ max = 1 corresponds 
to PM nom = rad (180°), we have 

PM nom = |0fast(7*Bcg ,nom 7 ^max ) | — £max ^ (29) 

While this is highly desirable, it is difficult to analyze the 
stability of the delay-perturbed system. This problem can 


be alleviated by using the Pade approximation to the delay 
transfer function, thereby obtaining approximately-linear- 
phase PM-gauges. In fact, the zero-error PM-gauge in the 
preceding section is the first-order Pade approximation to 
the delay function, but its phase function is quite nonlinear. 
Now, consider the third-order Pade approximation to (27) 


f-'fast^, £ ) — 


120 - 60T 0 £s + 12(7bes) 2 - (T 0 £s) 3 
120 + 60 T 0 £s + \2(Tq£s) 2 + (roes) 3 


(30) 


With To = u / C0 C g 5 nom and by Taylor expansion at e max = 0, 
we have 


P M nom — 1 0fast 0 ^cg ,nom 5 £max ) | 

(pmax^)^ 60(£ max ;r) 

120— 12(e max ^) 2 


2arctan 


- £n 


(31) 


Note that the estimation error of PM aom using the exact phase 
function |0f as tO&>cg,nom,£max)| is zero. The error of the linear 
approximation PM nom ~ e max n is given by 


^rr(£max) — 


2 arctan 


’ ( Tnax 71) 60(£ max 71 ) 

, 120-12(e max ^) 2 


■-1 


x 100% 


(32) 


which is a monotonically increasing function of e max maxi- 
mizing at e max = 1 to 0.64%, as shown in Fig. 5. 



Fig. 5. Approximation Error for the 3rd Order PM-gauge (30) 

Since practically PM nom is seldom more than -n rad 
(60°), for which the estimation error is given by 

err(£ max ) < err(^) = 0.00125% (33) 

which is accurate enough for any practical purposes. The 
test results of the approximately-linear-phase PM-gauge on 


the two nominal systems are presented in TABLE V. The 
approximation errors shown in the TABLE V are due to 
two sources: the approximately linear phase and the finite 
precision in implementing the transfer function (30), where 
by (32) and (33). The estimation error by the linear approx- 
imation is indiscernible, and most of the errors are due to 
the precision of the implementation. The SPM £ max values 
were obtained by solving (19) using the numerical method 
presented above using MATLAB. 

TABLE V 

Approximately-Linear-Phase PM-Gauoe Test Results 


Parameters 

Efast(-SS £) 

Case 1 

120- 108.8249ex+ 12(1. 8137ei) 2 - (1.81376s) 3 
120+ 108.8249M+ 12(1. 8137e.s) 2 + (1.8137es) 3 

Case 2 

120 — 68.99547o£s+ 12(1.1499es) 2 - (1.1499&S) 3 

120 + 68.9954r 0 ex + 12(1 . 1499es) 2 + (1 . 1499es) 3 


Param. 

SPMe m[lx 

£max ^ 

Exact PMnom 

£max7£ Err. 

Case 1 

0.3304 

1.0380 

71/3 = 1.0472 

0.8785% 

Case 2 

0.1660 

0.5216 

7T/6 = 0.5236 

0.3909% 


The third-order approximately-linear-phase PM-gauge 
(30) does increase the difficulty in evaluating the SPM £ max 
comparing to the previous two PM-gauge designs, but it 
is insignificant, and in fact worthwhile comparing to the 
simplicity and accuracy it provides in converting between 
SPM and PM. It is also noted that £ max T) = EmaxTt/tOcg.nom 
indicates how much transport delay the nominal system 
can accommodates before onset of instability, which is 
commonly known as the delay margin, which has been 
proposed in the literature as a candidate for stability metrics 
for NL systems. However, a practical method for assessing 
the exact delay margin has been lacking. The method 
developed here provides an analytical means to evaluate the 
(almost) exact delay margin for a given nominal system. 


IV. SUMMARY AND CONCLUSIONS 
A. Summary 

The main results of this paper can be summarized as 
follows. 

1) : Theorem 1 establishes a bijective relationship 

between the Singular Perturbation Margin (SPM) £ max and 
the phase margin (PM) of the nominal system with an 
approximation error on the order of £ 2 for monotonic-lag 
fast systems. In particular, the error will be zero if the gain 
of the fast system is identically unity, independent of £. 

2) : Lemma 1 reveals the relationships between the 
singular perturbation parameter £, PM of the perturbed 
system, SPM and PM of the nominal system, and the 
phase of the fast system at gain crossover frequency of the 
nominal system, which is an important result in its own right. 


3) : As an application of Theorem 1, a singular 

perturbation approach for time-domain PM assessment is 
proposed and illustrated by examples. The examples include 
design of PM-gauges and determination of SPM in the time 
domain from the state equation of the perturbed system 
consisting of the nominal and the PM-gauge systems. 
The effectiveness of this method is demonstrated with a 
zero-error PM-gauge design. 

B. Remarks 

1 ) : Relationship between SPM and PM of nominal 
system 

Theorem 1 establishes a bijective relationship between 
the SPM and the PM of the nominal system, with an 
approximation error on the order of e 2 , using a mono tonic - 
lag fast system. This relation not only makes it possible 
to assess the phase margin of the nominal system in the 
time-domain for LTI systems directly using a singular 
perturbation approach, but most importantly it establishes 
a basis for using SPM as a theoretically based, practically 
measurable and backward compatible metric for stability 
margins of LTV systems and NLTV systems. 

2) : PM estimation error analysis 

The PM estimation error is (9(e 2 ). Thus the error 
increases as £ increases. When the nominal system has a 
small PM, as shown Example 1 where PM nom is about 20°, 
the PM-gauge offers a result with very small error. When 
the nominal system has a larger PM, the fast system will 
need a larger £ max for its phase to reach the phase margin; 
therefore, the error will increase, as demonstrated with 
example 2. The 0(e 2 ) PM estimation error is due to the 
non-unity gain Mf ast (ft) C g, n0 m) of the fast (PM-gauge) system 
at ft) C g,nom- which causes the ft) C g,pt for the perturbed system 
to deviate from ft) C g,nom- Thus the error can be reduced by 
using PM-gauge with unity gain for all ft) (allpass filter) 
that is independent of £. 

3) : Assessment of SPM 

The first part of Section III offers a method to determine 
exact SPM in the time-domain from the state equation for 
LTI systems. For LTV and NL systems, general methods 
for determining exact SPM is unlikely to exist; therefore 
effective estimation methods should be developed in order 
to use SPM as a practical stability metric. For LTV systems, 
a recently developed differential algebraic spectral theory 
[17] -[22] could be used to assess the SPM in a manner 
similar to the aforementioned Chang transformation and 
multi-time-scale decomposition for LTI systems. The LTV 
SPM methods, once developed, could be extended to 
NLTV systems by virtue of the equivalence of exponential 
stability of an equilibrium state of a NL system and its 
linearization thereof. For NL system (including LTV systems 
as a special class), well-known Lyapunov function based 
singular perturbation results (see, for instance. Theorems 
11.3 and 11.4 in [11]) can be used to obtain conservative 


estimate of SPM. In this latter approach, a recent result 

on spectral Lyapunov function [23] would prove to be useful. 
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